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ABSTRACT 

A  linear  stability  analysis  for  a  two- layer  Binard  problem  is  considered. 
The  equations  are  not  self-adjoint.  The  system  can  lose  stability  to  time- 
periodic  disturbances.  For  example,  it  is  shown  numerically  that  when  the 
viscosities  and  coefficients  of  cubical  expansion  of  the  fluids  are  different, 
a  Hopf  bifurcation  can  occur,  resulting  in  a  pair  of  travelling  waves  or  a 
standing  wave.  This  may  have  application  in  the  modelling  of  convection  in 
the  Earth's  mantle. 
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SIGNIFICANCE  AND  EXPLANATION 


Flows  Involving  two  incompressible  viscous  fluids  exhibit  nonuniqueness 
in  the  sense  that  many  interface  positions  are  allowed  when  their  densities 
are  equal.  Two-fluid  flows  also  have  quite  different  dynamical  features  from 
one-fluid  flows.  The  one-fluid  s£nard  problem  in  which  the  fluid,  lying 
between  parallel  horizontal  plates,  is  heated  from  below  has  a  static  solution 
for  which  a  linear  stability  analysis  yields  no  complex  eigenvalues.  In  this 
paper  we  show  that  when  two  fluids  are  involved,  the  arrangement  in  horizontal 
layers  can  have  complex  eigenvalues  at  criticality  and  therefore  can  sustain 
disturbances  which  are  oscillatory  in  time,  nils  may  have  application  to  the 
theory  of  convection  in  the  Earth's  mantle,  which  is  sometimes  based  on  the 
assumption  that  convection  takes  place  in  chemically  uniform  layers. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
susMry  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


OSCILLATORY  INSTABILITY  IN  A  TWO-FLUID  BENARD  PROBLEM 
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1.  Introduction 

The  flow  of  two  immiscible  fluids  often  exhibits  phenomena  which  are  without  parallel 
In  the  flow  of  one  fluid.  An  example  is  the  steady  shear  flow  of  two  fluids  with 
different  viscoeities  but  similar  densities.  Such  flows  are  described  in  Joseph,  Nguyen 
and  Beavers1.  In  this  paper,  we  consider  the  Binard  problem  with  two  fluids  lying  between 
infinite  parallel  plates,  heated  frees  below,  and  we  look  for  new  phenoetena. 

In  the  Binard  problem  for  one  fluid,  the  'exchange  of  stabilities'  holds  and  all  the 
eigenvalues  of  the  linearised  problem  are  real.  In  the  two-fluid  problem,  we  have  both  a 
real  and  a  complex  spectrum. 

Busse  noted  that  convection  in  a  two-fluid  Binard  problem  heated  from  below  can 
admit  solutions  where  the  fluids  lie  in  layers  as  well  as  solutions  in  which  there  are 
convection  cells  of  one  fluid  surrounded  by  streamlines  of  the  second  fluid.  We  examine 
the  linear  stability  of  the  arrangement  where  the  fluids  lie  in  two  layers  with  a  flat 
horisontal  interface.  Zeren  and  Reynolds3  considered  this  problem,  including  the  effect 
of  a  linear  temperature  gradient  on  the. surface  tension  (Marangoni  effect).  They  state 
that  they  do  not  know  if  there  are  purely  Imaginary  eigenvalues  at  criticality.  They  note 
that  Sterling  and  Scrlven4  found  purely  imaginary  eigenvaluee  in  the  problem  where  the 
upper  fluid  is  inviscid  and  the  convection  is  induced  by  surface  tension,  which  depends 
linearly  on  the  temperature  of  the  free  surface  (Marangoni  problem).  Zeren  and  Reynolds 
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chose  to  compute  neutral  stability  curves  corresponding  to  zero  eigenvalues.  He 
concentrate  on  the  B£nard  problem  without  the  Marangoni  effect  and  show  that  the  equations 
are  not  self-adjoint.  We  give  an  example  of  a  situation  when  the  marginal  eigenvalues  are 
a  purely  imaginary  conjugate  pair  of  multiplicity  2  (the  same  eigenvalues  appear  for 
negative  wavenumbers).  Marginal  eigenvalues  of  this  type  are  associated  with  Hopf 
bifurcations  from  the  motionless  state  to  either  a  pair  of  travelling  waves  or  a  standing 
wave  (Ruelle5).  According  to  Ruelle5,  both  the  travelling  and  standing  waves  are 
solutions  to  the  nonlinear  problem.  If  they  are  both  supercritical,  then  only  one  of  them 
can  be  stable/  otherwise,  they  are  both  unstable.  The  possibility  of  travelling  waves  on 
the  interface  of  immiscible  fluids  may  have  application  to  the  modelling  of  mantle 
2 

convection  (Busse  ). 


2.  Umw  Stability  Analysis 


We  consider  the  linear  stability  of  a  two-dimensional  ( x* ,z*  >  problem  when  the 

*  •  • 

bottom  fluid  (fluid  1)  occupies  a  layer  from  r  >0  to  z  -  t 1  and  the  top  fluid 

ft  *  *  • 

(fluid  2)  lies  between  z  -  t 1  and  z  «  t  .  Asterisks  denote  dimensional  variables* 

•  *  •  t  *  • 

The  plate  at  z  -  0  is  at  temperature  Tg  +  AT  ,  AT  >  0,  and  the  plate  at  z  «  i  is 

•  * 
at  temperature  Tg.  Fluid  i  (i  »  1,2)  has  a  coefficient  of  cubical  expansion 

thermal  diffusivity  r^,  thermal  conductivity  ki#  viscosity  u1>  kinematic  viscosity 

* 

v^  and  density  at  temperature  Tg.  We  define  a  Rayleigh  number 

*  •  *3 

R  -  gdj At  A  /(x  w  ),  a  Prandtl  number  Pr  -  v} /*. ^  and  a  surface  tension  paraateter 
e 

Tn  “  St  /(x^),  where  S  is  the  surface  tension,  all  based  on  fluid  1.  There  are  6 
dimensionless  ratios i 

*  *  •  • 

m  -  M.,/u2,  r  -  Pj/Pj,  y  -  f/Kj.  C  -  k^j,  8  -  ®1/a2  "  *,/*  •  Denote 


We  choose  the  following  dimensionless  variables  (without  asterisks)! 

see  •  *2  *  *  *  *  **2  2  • 

(x,z)  -  (x  ,Z  )/t  ,  t  -  *tt  /t  ,  u  -  u  i  T  -  T  /AT  ,  P  -  p  t  /p^  where  Jl 

is  the  velocity  <u*,w*) ,  p*  is  the  pressure  and  T*  is  the  temperature.  The 

unperturbed  temperature  is 

Tg  ♦  1  -  A^z  for  0  <  z  <  t 

T-  (1) 

Tg  ♦  Aj(1  -  z)  for  tt  <  z  <  1 

where  Aj  -  «  "1  VT  and  A-  -  CA  and  the  unperturbed  motion  is  static.  A  linear 

*1  *  ^*2  ' 

perturbation  proportional  to  exp(ot  ♦  lax)  is  superposed  on  the  velocity,  temperature 
and  interface  position. 

The  perturbation  6  to  the  temperature  satisfies 

o8  -  wA1  ■  V2e,  for  0  «  z  «  *  , 

1  2  (2) 

08  -  wA2  -  ^  rl,  for  t,  <  z  <  1  . 

We  use  the  Boussinesq  approximation  in  the  Nevier-Stokes  equations.  Hence,  the  density  in 
the  buoyancy  term  is  approximated  by 

Pt(l  -  aA(T*  -  Tg) ),  1  -  1,2  (3) 


■*3- 


—mm  ms  YMVwaAfJV -V  --Ar-Sr.  «  -V  Jh  A  AA.V,  >  -  '-i..  • ..  ~-»  *—  *-»  ' 


to  yield i 


ou  ■  -Vp  ♦  RPr  8*  ♦  Pr  V2u  for  0  <  z  <  l, 

Ou  “  »r7 p  0e  +  —  Pr  V2u  for  i,  <  z  <  1 

-  r  B  -z  a  -  i 


where  es  is  the  unit  vertical  vector.  Incompressibility  yields 

V.u  -  0  .  (5) 

The  boundary  conditions  aret  u  ”  0,  8  »  0  at  z  •  0, 1.  The  following  linearized 
Interface  conditions  (see  Zeren  and  Reynolds3  for  coeiplete  derivation)  hold  at 
s  «  (.1  denotes  • 1  -  .j.  Continuity  of  velocity,  shear  stress,  temperature  end 

heat  flux  are.  respectively, 

[w]  -  [3w/3zl  -  0 
(y(32w/3s2  +  a2w)J  -  0 

(6) 

10]  -  h|XJ 
[k3  8/3*1  -  0  . 

The  kineautic  free-surface  condition  1* 

w  ■  oh  (7) 

where  the  perturbed  free-eurface  position  is  s  ■  *  h(x,t)  and 

h  ■  h^  exptiax  ♦  ot)  .  (B) 

The  conservation  of  volwe  of  the  incompressible  fluids  implies  that  the  h(x,t)  has  a 

zero  mean  value  as  a  function  of  x.  This  is  automatic  if  a  *  0.  There  is  a  difference 

between  a  ■  0  and  a  ♦  0,  the  fonwr  is  disallowed.  The  balance  of  normal  stress  is 

-  33w_/3*3  -  33w./3*3  ♦  3«2<1  -  -)3w,/3z 
■  7  i  mi 

t  (9> 

,  <~  -  1 )  ,  .  .  . 

-h«a<*{-4 -  +  tjAjd  -  jf)}  -  a  .Tn)  -  ~  (j  -  1  )3w^/3s  . 


We  will  show  that  the  above  problem  is  not  self-adjoint.  Hence,  the  eigenvalues  need 
not  be  real.  Let  0  be  a  strip  of  width  one  wavelength  2*/a,  covering  0  <  *  <  1. 

Let  be  the  part  of  fl  in  fluid  1  and  be  in  fluid  2.  Let  u  and  8  be  the 


ccaplcx  conjugate*  of  the  adjoint*  of  u  and  9.  the  asterisks  here  denote  the  adjoint 
and  the  overbars  the  complex  conjugates.  Integration  by  parts  of 

t  -*  2 

J  u  .  (au  ♦  7p  -  Pr  7  u  -  RPr9e  ) 

°r  ' 

♦  /  u  .  (—  u  +  7p  -  —  Pr  72u  -  9e  )  (10) 

a  *  r  *  b  •  rp  -* 

2 

+  /  5*  (00  -  wA  -  720)  *  J  -  (K[O0  -  wA]  -  720) 

rv  n  *• 


yields 


/  u  .  (ou  -  Pr[V2u  +  7(7. u  )]  -  A?  e  )  -  )  p7.u 

n, . ‘  n  ' 

+  /  u  .  (-  u*  -  —  t72u*  +  7(7.5*)]  -  0*e  ) 

a2  '  r  “  B  '  ‘  c  "* 

♦  J  0(o5*  -  RPrv*  -  725*>  ♦  J  0(^£  5*  -  v*  -  -i  V2ff*>  -  b 

n  ♦  ***  *• 


where  B  consists  of  boundary  integrals  taken  at  s  *  <1  over  one  wavelength  in  x.  We 
give  the  expression  for  B  later.  The  above  integration  is  facilitated  by  expressing 
/  72u  .  u  as  /  u  .  7.(7u  +  (7u)T),  where  superscript  T  denotes  the  transpose, 

n  -  -  n  " 

taxing  advantage  of  7.u  -  0,  to  obtain,  for  example, 

J  72u  .  u  m  j  u  .  72u  +  (u  .  7)(7.u  ) 

”(  *  *  °l' 

2L  .*  _«  -* 

*  J  Is  ‘57  +  +  2*  3l  -  u(a7"  ♦  3T>  -  *  • 

X-0 

(««t1 ) 


Choosing  u,  9,  p  and  its  derivatives  to  vanish  in  the  neighbourhood  of  the  interface,  we 

find  7  .  u  >0  and  other  adjoint  equations.  Since  7  .  u  «  0,  the  coefficients  of 

«# 

in  (11)  do  not  vanish  but  are  the  gradients  of  a  function  we  denote  by  p  .  Hence, 


in  fluid  1 


-«  -•  ?-* 

09  -  RPrw  -  7  9  »  0 


0  -*  Pr  „2-*  yh  s* 

7  a  -  -  v  u  ■  j  «  ^  ■  -7P  . 


in  fluid  2 


T  s*  RPr  -•  1  „2-* 

_  0e  .  _  w  --79  -  0  . 

We  examine  B  to  find  the  adjoint  interface  conditions •  The  integration  is  over  one 
wavelength  in  x  at  s  ■  t  . 


rr*,3ui  .  “a  *u2  , 


-#  -*  .*  -% 

u2  3u2  3w  3u^  3",  2  3w2 

♦  r  {w  ♦  st*'  -  tti *  3r> +  sr  - ;  w2  Jr 


3w#  3w*  39  9*  30  0  3ff*  jf* 

^  1  ,  2  „  2  i  .  jr*  0  1  2  2  .  2  2  B  1  . 

2wiir  +  :*2¥ri  +  ffi3r"r3r  +  r3i-“Bi3rdx 


a>«  ae  W 

da#  3w  3w  w  3w 

J  -  (Pl  -  2Pr  ^  w,  -  2Pr(w,  ^  ^  5^) 


#  3u,  3vt  ,  3u  3w^  3u,  3w, 

♦  Mis  1<3T ♦  sr> ♦  »,<:  ♦  3r‘  -  (3T  ♦  3T,)1 


.  30,  3f*  0,  35* 

^3T-  I91  371'  r3TJdX 


Conditions  (€)  and  (7)  yield 


;#cv;v 


3u  3v 

{u  ]  -  (9  ]  -  +  |T-)|  -  0  .  (' 

.  •*  .  * 

We  can  add  J  u  .  ?p  to  (10)  which  introduces  J  Ip  w]dx  into  B  in  (11).  Hence, 

a  * 

_• 

3w  3v  3w 

,  -*  wi  _*  2Pr  2  -*  1 

B  -  /  -  wt(p,  -  2VT-)  +  w2(p2  -  —  — )  +  wt(Pl  -  2Pr  57-) 


.  ___  35*  39*  9,  39* 

2Pr  2.  __r0  _ 1  2  _ 2,  . 

-  Vp2  -  ~  3r‘  •  "riei  5r  “  r  3rldx 


Condition  (9)  can  be  written  as 


*W1  1  **2  2 

p2  -  pi  ♦  2pr<ir  - :  3 r>  *  hI"i  *  •  v  * 0 


.  1 


Being  (6)  and  - 


M 1  -  BPr[-s - —  ♦  -  l]  and  -  PrTn  . 


2d  • 

a  5  — —  =  ,  the  laat  tent  in  ( 15)  is 

a*2  ,  . 


iJfliew,  -  I*"l«2)  • 


We  use  (6),  (15)  and  (16)  to  obtain 


r  -»  2wi  2Pr 

B  2Pr  5T>  ■  ,p2-  T-3T)J 


3*.  *,  . 

t5  j[p,  *  2»  it  -  ®,  e,  fly] 


m  i;*,9Mi-  *’*2"  •  9,  lr  ♦  r  ir  * 


V./ 


We  choosa 


[w*]  -  0 


and  uaa  (7)  to  find 


3»!  —  3w~ 


- 2pr ir  +  -?ir)  -  ( 


Wa  choosa 


a5*  3§*  35* 

IAJ  3T  1  ■  9i(dr  -  c  5T)dx 


35*  « 

TT  *  c  5T  "  0 


and 


_* 

3w 


3w_ 


)  -  (w>,  -  ;>2) 


aquations  f 12)— ( 14)  and  (17)-(19)  are  tha  adjoint  aquations. 


V  V' 


3.  Numeric*!  Scheme 

He  use  (S)  to  eliminate  u  so  that  in  each  fluid,  we  have  the  heat  equation  and  one 


■amentum  equation,  linear  in  0: 


2  2 

Pr(L  w  -  a  R9 )  -  OLw, 


wA  y  +  L0  »  o0 


Pr (—  L2w  -  0)  -  OLw, 

Bl  p 

wAj  +  ~  L0  ■  00 


for  0  <  z  <  i  1  , 


for  <  z  <  1 


where  L  =  32/3z2 


He  change  the  variable  z  to  z,  in  fluid  i  defined  by 


2  2 

z1  «  ■j—  z  -  1  and  tj  «  -j-  (z  -  1)  +  1  eo  that  the  z^  range  over  [-1,1]  in  each 

fluid.  He  then  expand  w(z^)  and  9(z^)  in  powers  of  Chebyschev  polynomials  T^z^) 

(Orszag6)  for  m  "  giving  a  total  of  4K  +  4  unknown  coefficients,  Together 

with  the  free-aurface  variable  hQ,  there  are  4N  +  5  unknowns.  There  are  6  boundary 

conditions  and  7  interface  conditions.  The  tens  of  highest  differential  order  in  the 
4  4 

■omentum  equation  is  3  w/3z  .  Since  we  choose  w  to  be  an  Nth  degree  polynomial,  the 

term  34w/3z*  is  of  degree  N  -  4  and  therefore  the  momentum  equation  is  truncated  at 

the  N  -  4th  degree,  yielding  N  -  3  equations  in  each  fluid.  Similarly,  since  the  term 

2  2 

of  highest  differential  order  in  the  heat  equation  is  3  w/3z  ,  we  truncate  this  equation 

at  the  N  -  2th  degree,  yielding  N  -  1  equations  in  each  fluid.  The  eigenvalues  of  the 

resulting  4N  +  5  square  matrix  equation  were  computed  in  complex  double  precision  on  a 
VAX 11-780  using  the  1HSL  routine  EIGZC. 

To  check  the  accuracy  and  convergence  of  our  computer  code,  we  computed  the 
eigenvalues  for  the  Benard  problem  in  one  fluid  with  Pr  ■  1,  R  -  2177.41  and  47005.6, 
a  -  2.  lt»e  eigenvalues  for  this  problem  are  real  and  are  given  by  Reid  and  Harris7.  The 

eigenvalues  at  criticality  (at  which  the  real  part  of  a  should  vanish;  are  less  than 

10~®  when  1  ■  15.  A  convergence  test  with  N  ■  15  and  20  showed  that  several  other 
eigenvalues  had  converged  to  at  least  5  figures  at  N  «  15. 


•„*  *  *  *,*  ‘  *,  *»  *  ,*•  .*•  .  '•  ,  ”  «  *  «  *  %  *  •  *  •  ‘  •  .  «*  «*  -•  .*  ,*  «*  /"  .  '  V  e_  <•  »  -w 
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The  computations  for  two  fluids  were  checked  against  Zeren  and  Reynolds'1  by  adding  an 

extra  tern  into  the  shear  stress  balance  at  the  interface  in  order  to  take  into  account 

the  Marangoni  effect.  He  define  a  Marangoni  number  based  on  fluid  1: 

Ha  «  (-  — -  and  our  shear  stress  condition  at  z  •  f.  is  modified  to: 

v  dT'  u.^  t 

2  2  2  2  2  2 
a  (m  -  1)w^  m3  Wj/3z  -  3  w^/3z  +  Ma.m.a  (0^  -  A^h)  »  0  . 

He  used  their  Table  2  for  the  values  of  the  physical  variables  at  16*C  for  benzene  lying 

above  water.  He  checked  our  eigenvalues  against  their  Table  3  for  t  ^  -  0.1  and  0.6 

for  heating  from  below.  Note  that  our  definition  of  the  R  and  Ma  are  different  from 

theirs.  At  -  0.1,  converting  their  parameters  to  ours,  they  find  criticality  at 

*  * 

Ma  -  1255.71,  R  -  178.3045,  c^AT  ■  0.00032537,  Pr  -  8.1,  a  -  3.5  and  Tn  -  460320.  He 
computed  (<J/Pr )  »  0.006186  using  both  N  «  15  and  20.  This  yields  0.00175  for  the 

eigenvalue  q  of  Zeren  and  Reynolds .  At  t  ■  0. 6,  their  parameters  in  Table  3  become 

*  * 

Ma  <*  4016.7153,  R  »  570.3736,  a^AT  -  0.0010408  and  0  -  2.5.  He  computed  -.00436  for 
their  eigenvalue  q  at  N  -  15  and  20.  In  both  cases,  we  also  found  stable  complex 
conjugate  pairs  in  the  spectrum. 
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4.  Numerical  «mbU» 

To  aid  tha  reader  in  the  interpretation  of  the  numerical  results,  we  recall  so<ss 

results  from  the  Benard  problesi  with  one  fluid.  In  the  sisiplest  case,  the  layer  is 

0 

bounded  at  *  -  0  and  z  -  1  by  atrasa-fraa  conducting  boundaries  and 


1  222f1  222 

-  (1  ♦  h)(n  «  ♦  «  )  ±  ij  (Pr  -  1)  (n  »  ♦ 


a2)2* 


a  RPr 

2  2  2; 
ns  ♦  o 


(20) 


for  n  ■  1/2/ see  e 

Hance,  for  a  -  0# 

o  -  -  j  (1  ♦  Pr)n2«2  +  I  |i»r  -  l|nV  <  0  ,  (21) 

end  as  a  ♦  «, 

2 

a  ~  -  f-  (d  ♦  Pr)  i  |Pr  -  l|)  <  0  .  (22) 

In  the  critical  case  o  •  0,  the  leaat  value  of  R  occurs  when  o  ■  ~  and 

2  2  3  2  /2 

R  -  (v  +  a  )  /a  .  These  f onaulas  are  for  stress-free  aurfacee  but  they  give  an  idea  of 

2 

the  variation  of  o(R,o  )  in  the  claeeical  case  of  one  fluid  between  rigid  boundaries. 

Now  we  consider  the  case  when  there  are  two  fluids  with  equal  properties.  This  would 
at  first  thought  appear  to  be  a  one-fluid  problem.  However,  it  is  easy  to  see  that  there 
is  a  solution  with  (0]  -  (a]  •  0,  hQ  i*  0  and  o  -  0.  He  shall  use  the  noaenclature 
introduced  by  Yih9  in  a  related  problesi  and  call  this  mode,  which  is  isiportant  when  the 
properties  of  the  two  fluids  are  different,  an  'interfacial  mode' .  He  track  eigenvalues 
as  we  vary  parameters.  Besides  the  interfacial  eigenfunctions,  we  have  other 
eigenfunctions  which  we  shall  call  Benard  modes. 

In  tracking  the  eigenvalues,  ws  shall  fix  all  the  parameters  so  that  there  is  a 
critical  <»  such  that  Be  o(a,R,Pr ,Tn,m,r,y,<; ,8,1^)  ■  0  with  Re  c  <  0  for  other  a. 

He  shall  exhibit  parameters  for  which  Xm  o  ?  0  at  criticality.  Hence  we  obtain 
oscillations  in  the  linear  problem  at  criticality  ('exchange  of  stabilities'  does  not 
hold)  and  the  nonlinear  problem  for  Benard  convection  in  two  fluids  can  have  time-periodic 


solutions  near  criticality 


Let  the  two  fluids  have  equal  densities  at  temperature  Tq  and  the  same  thermal 

diffusivities  and  conductivities:  r  •  Y  “  C  “  1.  Me  let  R  -  1695.7,  Pr  -  1 , 

*  • 

a^AT  »  0.001,  Tn  *»  0,  Ma  «  0,  m  -  1.1  and  6  -  0.9.  Thus,  if  fluid  1  occupies  the 
entire  flow,  the  Rayleigh  number  is  lower  than  the  critical  one  1708  (see  Reid  and 
Harris7).  If  fluid  2  occupies  the  entire  flow,  R  is  2072.52  and  the  flow  is  linearly 
unstable  for  a  range  of  a.  He  choose  i ^  »  0.4.  Figure  1  is  a  graph  of  the  growth 
rate  Re  o  against  a. 

He  are  approximately  at  criticality  when  a  ■  3.1.  In  this  case,  we  compute 
a  -  0.000072  t  15.9259  with  H  -  15,  20. 

The  five  numbers  next  to  the  curves  in  Figure  1  denote  branches  which  display 
different  features.  The  interfacial  node  ie  associated  with  branches  1,  3  and  5.  Branch 
1  can  be  obtained  from  the  interfacial  node  with  a  *  0  when  the  properties  of  the  fluids 
are  equal  (B  »  n  -  1)  by  moving  6  to  0.9  and  m  to  1.1.  nils  branch  is  real¬ 
valued.  Branch  2  is  associated  with  the  least  stable  of  the  Binard  nodes  for  a  single 
fluid  when  a  ♦  0.  This  branch  is  approximately  -9.87  at  a  ♦  o  and  would  correspond 
to  the  largest  value  of  (21).  Branch  2  is  real-valued.  Branches  1  and  2  coalesce  and 
split  into  conjugate  pairs  at  «  *  1.275.  At  a  -  6.79,  the  conjugate  pair  again  splits 
into  the  two  real-valued  branches  4  and  5.  Branch  4  is  associated  with  a  B&nard  mode  and 
remains  real,  decreasing  rapidly  as  a  is  Increased,  as  in  the  single  fluid  problem  (see 
(22)). 

Branch  5  is  an  interfacial  mode.  It  is  real-valued  and  negative.  The  stability  for 
large  a  which  is  associated  with  branches  4  and  5  is  explained  by  our  choice  of  6  and 
the  Boussinesq  approximation  (3).  He  consider  the  densities  p^l  -  a^(T  -  T0))  at  the 
unperturbed  interface  r  -  1 1  when  the  temperature  T  -  T0  is  given  by  (1).  Then  with 

A 

r  “  1  and  6  “  0.9,  we  find  that  Pj( 1  -  0.60^)  is  the  density  of  fluid  2  at  1 ^  -  0.4 

A 

and  P2<1  -  0. 54a 2)  is  the  density  of  fluid  1.  Hence  the  heavy  fluid  is  below  and 
gravity  nay  be  expected  to  stabilise  short  (large  a)  waves,  flie  interfacial  eigenvalue 
on  branch  5  is  discussed  in  {5. 


(VH0IS)TV3H 


F igur*  2  shows  ths  growth  rstss  versus  a  when  the  Rayleigh  number  is  increased  to 
R  -  2177.41  while  other  paraswters  are  fixed  as  in  Figure  1  and  shows  that  the 
instability  on  branch  3  is  associated  with  the  complex  conjugate  pair  of  branch  3  in 
Figure  1. 
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>totlc  teUyilj  of  the  Interfacial  Eigenvalue  for  Short  Waves 
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Vto  now  consider  disturbs ness  of  rapid  variation  whose  length  scale  of  variation  is  of 

th«  um  order  0(  1/a)  as  the  short  perturbation  wavelength10. 

Me  rescale  s  to  n  ■  a(z  -  t^)  and  let  n  be  0(1).  The  equations  In  fluid  1  are 

08  -  wA1  •  a2L  8  where  L*  5  32/3n2  -  1  and  <o  -  Pra2L  )L  w  «  -PPr8.  In  fluid  2, 

1  2  *  r  2  •  •  RPr 

08  -  w*2  -  —  a  L  8  and  (o  -  —  Pra  L  )L  w  »  -  — —  8.  The  Interface  conditions  are« 

w,  ■  w2  »  Oh  ,  (23) 

ISw^SnJ  -  0  ,  (24) 

{81  -[Alh  ,  (25) 

sJ^^/Jil2  -  S^j/Sh2  +  (at  -  1)w^  *  0  ,  (26) 


a2(^  33w2/3n3  -  a^/dn3)  +  3w1/3h.3a2(1  -  “>  ♦  ha3Tn 


-  hRa[-*—  +  SV1  "  rt>)  “  “  11  fc  5T 


Since  the  norsul  stress  condition  (2S)  contains  both  odd  and  even  powers  of  a,  ell  the 
variables  ere  formally  expended  in  powers  of  1/a.  2b  the  0th  end  1st  orders,  l  8  ■  0 
and  I»*w  •  0  In  each  fluid.  Using  conditions  (23),  (24)  and  (26),  we  obtain 
w1  -  C0(1  -  n)en  ♦  0(1/a)  and  w2  ■  C0(1  ♦  n)e”n  ♦  0(1/a)  as  a  ♦  •,  which  yields  to 
this  and  the  next  order,  3w/3n  *0  at  the  interface.  Hence,  the  normal  stress  condition 


(-  -  1 ) 

*(i  ^V*"*  ‘  »3'»,/8t’3)  ~  *  *[-5 - ♦  »2A2(1  "  rt>]  “  0 


where,  for  the  eoswnt,  surface  tension  has  been  neglected.  Tto  avoid  the  trivial  solution, 

2 

we  choose  0  «  oQ/a  +  0(1/a  )  for  large  a.  The  normal  stress  condition  yields 

.  (J-  «) 


.  _B —  (-E - ♦  i  a  (1  -  i->) 

0  2(-  ♦  1)  a, AT  22  ” 


r,.vfsJsivi^>.5S- . -v-V- v- Vv- 
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In  the  computations  (or  Figure  1,  the  asymptotic  formula  la  accurate  to  It  for 


a  *  20.  Me  confuted  -1.494  for  a  •  20  using  N  -  15  whereas  the  asymptotic  formula 
yields  -1.48. 

turning  now  to  a  consideration  of  the  effact  of  surface  tension,  we  find  that  when 
a2tn/R  -  0< 1 ) ,  then 


*0  * 


2<:  + 


i> 


(i 


i) 


*2*2 


( 29) 


We  computed  the  eigenvalue  for  the  peraawters  of  Figure  1  at  T  -  1  and  o  -  20  to  be 
-6.86  and  the  asymptotic  formula  yielda  -6.72.  Equation  (29)  shows  that  surface  tension 
ia  always  stabilising  for  short  wave  disturbances.  The  stabilisation  of  short-waves  by 
surface  tension,  even  with  adverse  density  ratios,  has  been  found  in  other  flows  such  as 
steady  shear  flows  with  two  immiscible  fluids  of  different  viscosities10'11. 
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